DESCRIPTION OF THE PRIMARY RELAXATION IN 
SUPERCOOLED LIQUIDS THROUGH THE TIMESCALE 

STEEPNESS FUNCTION. 



Valery B. Kokshenev^, Pablo D. Borges^, and Neil S. Sullivan^ 

^ Departamento de Fisica, Universidade Federal de Minas Gerais, 
CP. 702, 30123-970, Belo Horizonte, Minas Gerais, Brazil 
"^Physics Department, University of Florida, PO Box 117300, Gainesville, 

FL 32611-8400 
(February 2, 2008) 



Abstract 



The primary relaxation in glass forming supercooled liquids (SCLs) 
above the glass transformation temperature Tg is discussed in terms of the 
first-order (steepness) and the second-order (curvature) temperature deriva- 
tives of the observed primary relaxation timescale. We report new insights 
into the problem of the domain of the Vogel-Fulcher-Tamman (VFT) equa- 
tion, raised by Stickel et al. (J. Chem. Phys. 102, 6251 (1995), ibid. 
104, 2043 (1996)) and discussed by Richert and Angell {ibid. 108, 9016 
(1998)). A new ergodic-cluster Gaussian statistical approach to the prob- 
lem is given based on Onsager's thermodynamic principle. The primary 
relaxation is described by the VFT equation below the crossover temper- 
ature Tc (known from mode coupling theory (MOT)), and above Tc by an 
extended (VFTE) equation obtained after accounting for cluster-size fluc- 
tuations. The timescale is parametrized by a finite number of observable 



parameters such as the steepness function (my) at Tg , the MCT slowing- 
down exponent 7c, and the VFT and VFTE strength indices Dg and Dc- 
The latter are defined at Tg and Tc, respectively, for the strongly and mod- 
erately SCL states, which show absolute thermodynamic instability at the 
same VFT temperature Tq, associated with the Kauzmann temperature. 
For both states the limiting cluster-size characteristics are derived from 
experiment. A thermodynamic-dynamic correspondence is established be- 
tween the dynamic VFT equation and the thermodynamic Adam and Gibbs 
model. As a by-product, self-consistent equations for the excess heat capac- 
ity, entropy, and for the glass-transition characteristic temperatures Tq, Tg, 
and Tc are deduced. The problem of the "irregular" SCLs, which are not 
consistent with the standard VFT equation, such as salol, ortho-terphenyl, 
and bis-methoxy-phenyl-cyclohexane, is also discussed. 
PAC Nos: OO.OO.+g, OO.OO.Cn, OO.OO.Rh OO.OO.+g, OO.OO.Cn, OO.OO.Rh 
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I. INTRODUCTION 



The understanding of noncrystalline liquids and the formation of glasses by a dynamical 
slowing- down, and a rapid increase in viscosity at the glass transition temperatures Tg (near 
2/3 of the melting point T^) is currently seen as a major intellectual challenge in condensed 
matter physics |T|. There is therefore intense interest in studies, experimental and theo- 
retical, including simulations, of the temperature dependence of the transport properties of 
supercooled liquids (SCLs). In particular, data on the primary relaxation times r-^^^^'* that 
are determined from dynamic experiments that study viscoelastic, dielectric, conductivity, 
light scattering and mechanical relaxation data of glass formers 0, offer one of the main 
keys to the problem 

Unlike the slow secondary relaxation process which has a constant activation energy 
^ 2AkBTg (established in Ref. for the Arrhenius relaxation, where ks is Boltzmann's 
factor), the thermally equilibrated SCL states, observed within the temperature domain 
Tg < T < Tm, exhibit a pseudo-Arrhenius form, i.e., tt = Too exp(i?r/A;BT), and are dis- 
cussed traditionally in terms of the SCL activation energy Et- The most widely employed 
Vogel-Fulcher-Tammann (VFT) timescale fitting phenomenological form with the activation 
energy Et = B /{T — Tq) is introduced by the fictive critical temperature T^r (= Tq, VFT 
temperature). Tq lies below the thermodynamic glass transition temperature Tg established 
by thermodynamic experiments such as scanning calorimetric or volume-expansion data. 
Analysis of the temperature-derivative behavior of the a-relaxation timescale tt given by 
Stickel et al. in Refs. [0-^, as improved and extended by Richert and Angell's analysis of 
the thermodynamic data in Ref. , determined the real temperature observation domain, 
or observation window, of the VFT form. Overall tt fitting forms, distinct from the VFT, 
were proposed within the whole SCL observation window by Kivelson et al. in Ref. 11 or by 
Schulz in Ref. fl^ that, respectively, admits the "avoided critical points" |1^ {Tcr > Tm), 
or rules out the non-zero critical temperatures {Tcr = 0). 

Within the context of first-and-second-order phase-transition approaches to the problem 
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of the liquid-to-glass transition (see Refs. |[T^-p3|) the dramatic increase of the relaxation 



times t!^^^^ observed in SCLs upon cooling, is due to the growth of underlying correlation- 
length scales. The lengthscale peculiarities, even when being experimentally established, 
do not often lead to a clear physical understanding. In the case of the VFT form the 
physical meaning of the experimentally determined VFT temperature Tq, is justified by its 
closeness to the thermodynamic Kauzmann temperature T^- At this critical temperature 
(Tk = Tcr) the excess liquid-over-solid SCL entropy, or configurational entropy, vanishes. 
This property is thought to be associated with the absolute thermodynamic instability of 
the thermally equilibrated states of the SCLs with respect to liquid-to-solid rearranging 
processes that extend from the ergodic region {Tg < T < T^) to lower temperatures. The 
theoretically motivated ergodic instability (caused by divergence of the zero moment of the 
relevant relaxation-time distribution function [jl9|) establishes the thermodynamic-dynamic 
(TD-D) correspondence equation, T^^^ = Tq^"* = Tq. This equation was experimentally 
justified by Angell |^ within a certain range of experimental scatter which often lies ||2^ 
inside the error of experimentally measured characteristic temperatures. 

Thermodynamic models for SCLs are commonly tested by the t^^"^^ fitting analysis. The 
VFT form has proved to be ineffective for establishing simultaneously all three adjustable 
parameters (too , B and Tq) within the overall fitting analysis of the observed SCL vis- 
cosity (see discussion in Ref. An improved analysis for the tested forms was proposed 
in Refs. on the basis of fitting experimental data r^'^^'^-'and additionally their first- 
order {dlogr!^^'^^ /dT ) and second-order {d"^ log t!^^"^^ / dT"^) derivatives. These temperature- 
derivative analyses permit one to exclude the unknown fitting parameters, step-by-step, by 
making the three-parameter fitting procedure self-consistent. As a result, the validity of the 
VFT form, when tested within the whole ergodic relaxation range, was questioned for all 
SCLs investigated above the so-called a-/3-relaxation bifurcation temperature Tb- Notably, 
this crossover temperature was often found close to other known crossover temperatures 
and Tc introduced, respectively, as a scaling temperature by Rossler and co-workers P3 



and as the dynamic critical temperature of the mode coupling theory (MCT). Within these 



studies one can admit that to a good approximation, Tb = = Tc (see e.g. Table II in Ref. 
0)- 

Physically, the MCT crossover temperature Tc signals structural transformations, or 
liquid-solid-like restructuring in SCLs, from moderately supercooled states of "weakly cou- 
pled fluids" (Tc < T < Tm) to strongly supercooled states of "strongly coupled fluids" 
(Tg < T < Tc) as discussed in Ref. The analyses given in Refs. [f7|-[lOt| in terms of the 
first order derivatives (related to the timescale steepness) and the second-order derivatives 
(related to the timescale curvature) of the temperature dependence, clearly distinguish the 
two regions, similar to the case of the MCT theory [4], below and above the dynamical 
crossover temperature Tg, associated here with T^. The data for moderately SCL (or weak- 
SCL) states in all glass forming liquids that have been studied are in a good quantitative 
agreement with the modified VFT form given as the high- T VFT (HT- VFT) version in the 
old canonical VFT form but with new adjustable parameters, including the HT-VFT char- 
acteristic temperature Tq^'^\'^ Tg, see Table II in Ref. [§). This fundamental finding is 
exemplified by the case of the simple molecular liquid propylene carbonate (PC) in Fig.l. As 
seen from the analysis given in Fig.l, the validity of the canonical VFT form is limited by the 
temperature domain Tg < T < Tc, corresponding to the strongly SCL states, characterized 
by relaxation times 10^'^ < r^^^^ / s < 10^, or by viscosities 10^ < t]!^^^^ /P < 10^^. 

Another interesting description of the SL states is also provided by analysis of the ther- 
modynamic (heat capacity) data that was elucidated by Richert and Angell [jlO[ on the basis 



of the configurational entropy Sc data. This data is schematically plotted in Fig. 2. 

Careful complex dynamic-thermodynamic analyses enabled one to establish a physical 
meaning for the standard VFT fitting form, that was justified in terms of the thermodynamic 



Adam-Gibbs model (ACM) introduced in Ref. |T^. The VFT form is therefore treated as the 



VFT equation that provides a good description of the observed dynamic and thermodynamic 
relaxation of the strong-SCL states (see AG- line in Fig. 2). Within the context of percolation 
theory [Q, Tq may be associated with the critical temperature where the correlation length 
of the strong-SCL states diverges in accordance with the thermodynamic instability at the 



Kauzmann temperature T^. A microscopic version of this kind of instability was given 
by Bender and Shlesinger Ref. by invoking a defect-aggregation mechanism that was 
recently justified [|28[ for PC by the pressure-dependent experiments. 

In reality, strong-SCL ergodic states do not survive below Tg (see solid squares in Fig.2). 
As shown in Ref. I^Uj, this is due to the ergodic-to-nonergodic crossover from the short- 



range strong liquid correlations, characteristic of solids, to somewhat long-range weak liquid 
correlations, characteristic of gases. Conversely, no physical justification was given for the 
HT-VFT form introduced in Refs. [jTHToH. The main aim of the current study is to provide a 
physical description of the weak-SCL states on the basis of both the extended thermodynamic 
AGM and the dynamic VFT equations. 

The old idea that anomalous temporal (non-Debye) and anomalous temperature (non- 
Arrhenius) behaviors related to the SCL relaxation, may be attributed to rearranging regions 
or domains ||21|, or droplets ||20| or clusters p9[ of different size and structure, is still 



widely explored and remains fruitful. There is today considerable direct evidence for the 
existence of these kinds of regions. They have been observed on a microscopic scale by 
multidimensional NMR, photobleaching, excess light scattering and dielectric hole burning, 
as reviewed in Ref. |T]], and by X-ray scattering data, discussed in Ref. The percolation 
theory, improved by mesoscopic-domain energy-size scaling, as shown by Chamberlin in Ref. 
113011 , provided a description in detail of the observed relaxation function. An alternative 
percolation-theory approach, combined with the relaxation-time cluster-size scaling law, 
proposed in Ref. shed light on some universal features of the primary relaxation near Tg 



in the structural and orientational glasses [Q. Additional cluster-approach studies of the 



dynamical spectra with accounting of the MCT scaling forms, explained the mechanisms 
for the primary relaxation in polymers and those of the secondary relaxation in SCLs. 
Furthermore, exploration of Stauffer's cluster-size scaling form, based on the percolation 
theory, made it possible [^] to justify the underlying ergodic- nonergodic transition that 
emerges at the ergodic crossover temperature Te (< Tg), where the cluster-size Gaussian 
distribution of the strong-SCL states is transformed upon cooling into the Poisson form (see 



also discussions in Refs. |1T9| , |35| ). 

In the current study we extend the cluster approach to the problem of relaxation of 
ergodic states in SCLs. We will demonstrate that the structural heterogenous fluctuations, 
being small in the case of the strong-SCL state, are responsible for moderately SCL behavior 
as revealed in Refs. [|7|-pi][] by deviations from the standard VFT fitting form. The paper is 
organized as follows. In Sec. II we discuss phenomenological descriptions of the observed 
primary relaxation timescale in terms of the known fitting forms and macroscopic parameters 
commonly used by experimentalists. Instead of the frequency observation windows, given 
by 1/t!j^^^\ which are distinct for different dynamic experiments, the timescale is introduced 
through the temperature domains for the SCL excitations established with the help of the 
characteristic temperatures. In Sec. Ill we give the dynamic and thermodynamic description 
of SCLs based on the statistical approach to the rearranging regions introduced by Adam 
and Gibbs. By accounting for their size fluctuations, we develop extended versions for the 
AGM and the VFT equations. The qualitative and quantitative analyses of the proposed 
physically justified fitting forms are given in Sec. IV. The conclusions are summarized in 
Sec. V. 

II. MACROSCOPIC PARAMETRIZATION OR THE PRIMARY RELAXATION 

TIMESCALE 

A. Strongly Supercooled Liquids 



It was experimentally proved by Richert and Angell in Ref. 0] that the widely employed 
AGM equation 

log.rr'^A+^ (1) 
and the phenomenological VFT form given in the two equivalent forms, namely 

loe T^^^^^ - loe r(^^^) + ^ (2^ 



rr"' = rr"'expj^ (3) 

are self-consistent when subject to the conditions A = logr^^^^ and Tq = Tk ■ These 
equations apparently work well in describing the strong-SCL state relaxation in terms of the 
MCT within the temperature domain Tg < T < T^. This domain was studied [|1^] simul- 
taneously on the basis of dynamic and thermodynamic experimental data, more precisely, 
through the primary relaxation timescale T!f^^\ deduced from the dielectric loss spectra, 
and the configurational entropy: 

'Tk 



S,{T) = ASt = I ACrT-^dT (4) 



with the excess liquid-over-solid isobaric specific heat capacity ACt = Cp^^'^^ — Cp^°''\ Ac- 
counting for the fact that at high enough temperatures |26|,0 AC^^^"* oc T~^, the configu- 



rational entropy (HI) was given [|T0|j3H] in explicit form, namely 



S,{T) = ASoo (l - ^) with AS^ = ^ (5) 

that establishes useful relations between the observed model parameters. It has been demon- 
strated [|l^ that unlike the VFT form, restricted by the strong-SCL states, Eq.(|^) does a 
good job within the whole range of the thermally equilibrated states (see solid line in Fig.2). 

The theoretical foundation of the experimental findings by Richert and Angell are not 
limited by the aforesaid TD-D correspondence equation [|T0|, i.e., T^^^ = Tq^^ = Tq. 
One can see that the existence of the strongly SCL states, based on the dynamic VFT 
equation and the thermodynamic AGM within observation windows, evidently requires in 
addition Tg^^"* = Tg^^ = Tg. Here the characteristic temperature Tg^^^ is introduced as 
a thermodynamic parameter derived from the maximum of '"^(Tf ^^), observed by the 
scanning calorimetry, and the corresponding laboratory glass-point Tg^^ that is deduced from 
the dynamic experiments and defined by the equation r^'^^^-* = lO^s given at T = T^^\ 
Unlike the case of the first TD-D equation, there is no systematic study to verify the precision 
of the second equation. Except for the case of SCLs (see e.g. Fig.2 in Ref. ^9|) there is even 
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no accepted international convention on the definition of tlie tliermodynamic value Tg^^^ 
for the various glass forming materials. 

Furthermore, Eqs.(|l], ^ cannot be extrapolated to high temperatures beyond their tem- 



perature domain. Hence, the justification given in Ref. |10] for the VFT and the ADM 
corresponding macroscopic equations, i.e., B = C/ASoo and A = log Too , should be re- 
considered. The same equations might be introduced by application of the ergodic hypothesis 
employed for SCLs in Ref. This means that the strongly SCL states, that emerge above 
the crossover-ergodic-temperature (see Eq.(7) in Ref. [^]), can be introduced through 
certain observables estimated by thermodynamic and dynamic averages. Physically, this cor- 
responds to the fact that the thermally equilibrated states are equivalent to ergodic states 
observed in the dynamical experiments. If this is true, this requires new equations for the 
ergodic time scale and length scale given, respectively, by 

In other words, Eq.(^ introduces the minimum magnitudes for temporal and size charac- 
teristics necessary employed for the mesoscopic description of the SCL states. 



B. Moderately Supercooled Liquids 

In contrast to the second-order thermodynamic transition approach, MCT treats the 
liquid-to-glass structural transformation as a dynamic process driven by liquid density fluc- 
tuations. In the idealized version of MCT [4] the ergodic-to-nonergodic transition occurs at 
a unique crossover temperature Tc(= T^r) at which neutron spin scattering measurements 
signal the appearance of the glassy-like correlations given by the so-called non-ergodic pa- 
rameter, which is analogues to the Edwards-Anderson spin-glass order parameter in the 
case of weak-SCL states. Moreover, the MCT predicts |^ the algebraic-type scaling for the 
primary relaxation timescale, namely 



Tr 



(MCT) 

{MCT) _ Tqo 



(T - Te 



for T,<Ty<T <Ta<T^ (7) 



with the slowing- down exponent 7c. This observable is rather sensitive to the location of 
Tc as well to the temperature region of the fitting experimental curve P,|7|,p|, p!T| , ^ . The 



temperature domain of Eq.(|^ is much narrower than that for the case of the HT-VFT 
form given by < T < Ta, where Ta is the Arrhenius crossover temperature. A 
temperature Ty is introduced in Eq.(|^) to designate the lowest temperature point where the 
algebraic MCT form fits experimental data (see Fig.l). Unlike the first estimates reported 
for the slowing-down exponents that exceeded magnitudes of 5.5 in SCLs (see Tablel in Ref. 
pop, the values established by the MCT Eq.(|^) in glass-forming materials are in the range, 
2 < 7c < 3 for liquids Eq.(^, and 3 < 7c < 4 for polymers. For macroscopic observable 
parameters of the timescale in glass forming molecular liquids and polymers see, respectively, 
cases A and B in Table 1. 

C. Fragility and Characteristic Temperatures 

One of the salient features of the primary relaxation is the so-called time-temperature 
superposition principle [^]. Its physical justification has not yet been established but it 
is employed effectively to scale the dielectric loss spectra Meanwhile, analyses of a 
number of the a-timescale characteristics near Tg for several experimental studies of glass 
forming materials |@], have revealed a number of underlying tempera/- sj>atm/ correlations in 
explicit form. The analysis was given in terms of the phenomenological relation, i.e., 
nig = (250 — 320)Pg expressed through the fragility nig that is related to the steepness 



index [|l0|, and the Kohlrausch-Wilhams- Watts (KWW) stretched exponent (3g. 

Incorporation |3^ of the growth-correlation scaling law, employed in the order- disorder 
kinetic transition theory into the relaxation function of the SCLs led to a description 



of the aforesaid temporal-spatial correlations for the dynamical exponents in a form 
(3g = {1+ rrig/mQ)^^ where is a universal parameter (because it is material-independent) 
within a certain glass family. In addition, the equation for the growth-correlation dynamic 
exponent, i.e., Zg = mg/SniQ , was established in Ref. [33] (see also proposal by Ngai in Ref. 
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ii)- 

Within MCT, the characteristic temperatures, as well as the dynamical exponents, are 
treated as certain equilibrium thermodynamic parameters, which, as demonstrated in Ref. 
are not independent. This statement is exemplified by the widely employed equation 



^ = with ml = log.o i , (8) 

To nig-m* ^ Ta 

where m* = 16 ± 2 is the universal parameter known as the lower fragility limit 0. As 
follows from Eq.(|D, the uncertainty in the m*, caused p3|j4^ by deviation of the inclination 



of the straight lines in Fig. 3, is due to the uncertainty of the timescale unit parameter 
t!^^'^^ = Ta = 10~^^''^"'^s, given by Eqs.@, and to that of the glass-transition characteristic 
time Tg = T* = lO^^^s. Eq.(D was first introduced in Refs. on the basis of the VFT 

equation @, namely 

rf = rir^^) exp ^ with st = ^ - 1 for Tg < T < T„ (9) 

where the timescale auxiliary function Et is introduced. 
Let us define the timescale steepness function by 



with 



for the particular case of Eq.(|^). Eq. (pT]) , estimated at T = Tg, is known as the 
fragility that was proposed by Angell for the strong-to-fragile classification of glass 
forming materials We see that unlike the case of r^^^"^\ the function m^^'^^ has only 
two free parameters, Tq and Dg. Moreover, the latter can be excluded from the two-equation 
system given by rj^'^^-'and m^^^^\ using Dg = m*eglnlO, that eventually results in the 
known Eq.®. 
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As shown above, Eq.(P) is a formal consequence of the self-consistency between the VFT 
equation t'^^'^^ and its first derivative given by the observed fragility ■m^^^'^\ Physically, the 
existence of Eq.(|[) is due to the presence of the strong-SCL states at the lowest temperature 
Tg parametrized by these equations. One can therefore introduce a similar relation (see Eq.8 
in Ref. ||2^), namely 

^ = ^^^^^ with ml = log.o ^ (12) 
To nig - m* Ta 

associated with the weak-SCL states at the lowest temperature T^. The validity of Eq.(|12|) 

was established numerically with the universal parameter m* = 7± 1 (see Fig. 2 in Ref. p3[ ) 

that corresponds to the adopted value r* = lO^^^s. Eq.(|T^) was deduced by application 

of the Taylor series for the steepness function (|TD|) near Tg and Tc along with the following 

assumptions: (i) in contrast to the case of the timescale the steepness function rriT is 

analytical and smooth beyond Tq; and (ii) similar to the strong-SCL state, described by 

nirp with parameters Dg and Tq, the weak-SCL states can be described within the scope 

of the HT-VFT form with help of m^^^\ using D^^^^ = m* Sc In 10, and ec= Tc/Tg — 1, 

i.e. by the formal substitution of Tq and Tg in the standard VFT equation by Tg and Tc 

in the HT-VFT form, respectively. It is remarkable that the idea to introduce in Ref. |2^ 



the steepness function mrp on the basis of the HT-VFT analytical form was stimulated by 
the uncertainty of the fitting analyses of t^^"^^ discussed by Kivelson et al. in Ref. JIT]. 



Experimentally, the HT-VFT form was established by the first-order derivative analyses 
given by Stickel et al. in Refs. [0-0]. 

The most striking consequence of the existence of the two distinct SCL states described at 
Tg (in the vicinity but below Tc), respectively, by Eqs.(|D and (|12D is an inference on the self- 
consistency of the principal characteristic temperatures. Indeed, excluding the fragility nig 
from Eqs.(^ and (|12D one obtains the following characteristic-temperature equation (CTE), 
namely 



- c 



Tc Tg(^ + l)-T, 
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(13) 



(that was analyzed analytically in Ref. under the model simplification m* = 2m*). One 
can see that beyond the scope of this simplification, Eq.(|T3|) includes the particular case of 
the MCT ideal liquid-to-glass transition that occurs at a unique crossover temperature given 
by the equation Tg = Tq = Tc. 

We have shown that representation of the CTE in the form given in Eq.(^) is ensured 
by application of the VFT fitting forms to both the strong and the weak SCL states. On the 



other hand, there exist a number of fitting forms (see e.g. Refs. P7|j28|| and for a short review 
see Ref. 0) and therefore the question arises as to whether the fact of self-dependency of 
characteristic temperatures given through the CTE is also limited by exploration of the VFT 
forms. 

Using the Lindeman melting criterion for the mean-square displacements of pre-transition 



solid-like clusters in SCLs, Novikov et al. in Ref. deduced the CTE in the "square-root" 
form: Tg"^"* = \^TqTc. Meanwhile it was shown ||2^ that Tg°^^ and the "cubic-root" CTE, 
i.e., Tg"^"* = ■^a/TqT^^, are respectively, the lower and the upper limits of the solutions of 
Eq.(p!3D considered under the model constraint m* = 2m*. Physically, this means that the 
experimental data on Tg'^^^'* and theoretical predictions for T^, deduced from Eq.(^3|) under 
conditions Tq = T^"^^^ and Tc = Tjt'^^\ are expected to satisfy the equation Tg = Tg'^^^^ 
within the experimental error of ±5% (see the last column in Table 1 in Ref. ||23|) and, in 
addition, to satisfy the inequalities Tg'*""-* < Tg < Tg'^^\ 

Odagaki proposed a controlled-diffusion approach to the SCL glass-transition problem in 
Ref. |J19||. A glass transformation process was characterized by the critical VFT temperature 
To (= Tk) and the crossover temperature Tc (= T^) related to the peculiarities of the pseudo- 
divergent a-timescale, and introduced through the true divergence of the first moment 
of their distribution function. As a consequence the CTE in the form TcASc = 2TgASg was 
deduced. Taking into account that the excess entropy ASt, given in Eq.(^), in real SCLs is 
well approximated by Eq.(j|), we represent here the Odagaki CTE in the following equivalent 
form, namely 
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Again, one can see that the aforementioned MCT ideal transition is also satisfied by Eq.(|l^). 
A comparative analysis of the distinct CTE form is given in Fig. 4 for glass forming liquids 
(case A) and polymers (case B). As seen from Fig. 4 the Odagaki CTE prediction (p!4D is very 
close to that given by Eq.(p!3D with a value for m*/m* that lies within the experimentally 
established uncertainties, i.e., m* = 16 ± 2 and m* = 7 ± 1. 

III. DYNAMICS AND THERMODYNAMICS OF SUPERCOOLED LIQUID 

STATES 

A. Cluster Statistics 



Solid-like clusters were treated in Refs. p9| , ^p4| as short-range random-size correlated 



regions of SCLs that permits one, at a mesoscopic level of consideration, to employ the 
percolation [^,^ or controUed-diffusion p3| , ^ models. Within this context, clusters are 



equivalent to solid-like domains [p!3| , p0| , or droplets ||2^ or CRRs (cooperatively rearranging 



regions) as introduced in Ref. ||14|. A more detailed specification of clusters depends on 
ways of incorporation of short-range correlations in the description of their evolution. 

We treat clusters as conglomerates of z molecules that are characterized at a given 
temperature T and pressure p by the solid-state chemical potential n^"''^ (p) . These finite-size 
solid-like clusters are in thermal and mechanical contact with the remaining N — z molecules, 
which are in the normal liquid (NL) state and characterized by a chemical potential IJ''t''\p)- 
The typical cluster, or characteristic cluster, is in dynamic equilibrium with SCL, and is 
introduced by the mean-size configurational average, i.e., < z >c= Zt{p) through the cluster- 
size distribution function Pt{p,z). For the thermally equilibrated SCL states we admit 



the ergodicity hypothesis ||29|, i.e. the configurational averaged values are equal to those 
estimated at thermodynamic equilibrium for temperatures T > Tg > Te, above the ergodic- 
nonergodic transition temperature Te [29]. 
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The isobaric process, as a model mechanism of 2;-cluster formation, can be described 
by variation of the Gibbs (thermodynamic) potential 6G^^°''\T,p, z) = —S!^''''^6T + ij,ip°''^6z. 
Solidification, or cluster growth, if assumed at the thermal equilibrium temperature Tgg, is 
defined by minimization of the whole system potential, that gives rise to the requirements 
imposed on the chemical potentials, i.e., = /xig^'', and on the entropies, with S'eg°''' < 

Seg'^\ given at T = T^g (see e.g. Fig. 15 and discussion in Ref. |^). Stabilization of 



the liquid-solid boundary upon decrease of temperature below the equilibrium temperature 
is also ensured by minimization of the whole Gibbs potential and therefore by positivity 
of the excess chemical potential, i.e., AfiT = — /i^*'^'' > and the mesoscopic- scale 



configurational entropy, i.e., ASt{z) = s!^^'^\N — z) — s!^°''\z) > at T < T^g (hereafter 



we 



omit the pressure p as an irrelevant parameter). Note that the macroscopic configuration 
entropy introduced by Eq.(^ is defined here as ASt =< ASt{z) >c- The latter is also 
expected to be valid far from thermodynamic equilibrium. 

Within the fluctuation mechanism for the formation of heterogeneous clusters, the prob- 
ability for the formation of a solid cluster of size z, is formally given by the probability 
PT{z)dz that the thermodynamic variable z lies between z and z + dz. This probability is 
driven by the temperature-dependent whole-system entropy namely 



St_{z) 



Pt{z) = exp — ^ . / — ^ exp 



2Az\ 



(15) 



Here the characteristic cluster of size zt is in thermodynamic equilibrium with the thermal 
bath, and therefore the entropy of the whole system St{^z) = s!^"'''' + s!p'^^ + S"^™^-* (which 
includes the solid- liquid mixing term s!^^^'') has a maximum aX z = zt- The cluster- 
size fluctuations Az\ =< {z — ztY >c are introduced formally by the relation Az^ = 
— {d"^ St / ksd"^ z)~}^^^, and Eq. (p!5|) is expected to work near equilibrium, i.e., for T ^ Te,. It is 
notable that the real SCLs, studied in dynamic experiments, are not close to thermodynamic 
equilibrium, but the Gaussian size-distribution function Pt{z) of Eq. (p!5| ) was experimentally 
tested and justified by Chamberlin et al. in Refs. PD|,P5| for a number of SCLs. 



Let us treat the formation of a cluster of size z at fixed pressure p and temperature 
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T as an isobaric-isothermal fluctuation process. This means that sohd-hke clusters are not 
in thermodynamic equihbrium with the fluid, i.e., T < T^q, [T refers to the solid-state 
subsystem and T^q to the thermal bath) but they appear and disappear from a fluid volume 
of size z with a CRR frequency r^^{z), where tt{z) is the lifetime of a z-cluster. The 
probability of formation of this kind of cluster is given by TgT^(z) and is determined from 
Onsager's principle, i.e. by the aforesaid mesoscopic size-scale entropy ASt{z), namely 

W^^iniz)' 



Kz) = exp ( ) = exp 



ksTeq 



(16) 



Here IVmm = AGj'(z) corresponds to the minimum work (see e.g. Eqs.(20.8) and (114.1) 
of Ref. |^8|) required for solidification of z molecules driven by cluster-size fluctuations. 



Within the isobaric-isothermal mechanism, proposed by Adam and Gibbs [|IJ], one has 

As a result the exact result of Eq.(|l^) can be extended to 
non-equilibrium states with T ^ Teq through the relation 

Mz) = r^e.J^) (17) 



which introduces the cluster mesoscopic- scale relaxation time. This in turn is related to the 
observed macroscopic timescale, namely 

POO 

tt =< Mz) >c= / TT{z)PT{z)dz. (18) 

Within the proposed cluster statistical treatment, the corresponding macroscopic lengthscale 
is given by 



zt =< z >c= / zPT{z)dz. (19) 
Jo 

On one hand, these two scales, experimentally justified above the glass transition temper- 
ature, are ensured by the thermodynamic-dynamic correspondence principle. On the other 
hand, the cluster-formation mechanism, specified by Pt{z) through the observables given in 
Eqs.([T8|) and ([T9|), should satisfy the constraints that are naturally required by the macro- 
scopic and mesoscopic scales. In other words, the observables should obey the conditions 
Tt > rj™™-* = Ta and zjn > z!^™^ = Za where the limiting mesoscopic-scale magnitudes 
and Z(i circ estabhshed in Eq.(||). 
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B. Adam-Gibbs Model 



The Adam-Gibbs molecular-kinetic theory |n| was based on Eq.(]IBD given for the coop- 
erative transition probability for rearranging regions of z molecules with the minimum work 
IV^ff ^ = A/i(^^)2 that was estimated for the isobaric-isothermal ensemble of equivalent 
thermodynamic subsystems, where Afx^"^^^ is a solid-over-liquid excess chemical potential, 
approximated by a constant. The average transition probability follows from the 

AGM macroscopic timescale, namely 

rr^=rr^exp ( ""'Tf' ]- (20) 



The smallest GRR of a "critical size" z!^'^\ that still permits a transition, corresponds to 
the lowest critical configurational entropy for one GGR, s* = ks In 2, introduced as the first 
estimate in Ref. [|14|. The configurational entropy Sc{T) of the whole system given in Eqs. 



(I) and @ for one mole of the SGL (A^ = Na), consists of Na/z!^^^ of self-similar smallest 
GRRs and is therefore estimated here as 

= s:^. (21) 

One can see that Eq.(|2T|) is equivalent to Eq.(20) in Ref. [g. Eqs.(|2D|) and (^ deduced 
above result in the known Eq.(0) with the AGM parameters, namely 

^^,og...r'andC^^^^^^. (22) 

Kb 

Note that the AGM was motivated to validate the phenomenological Williams-Landel- Ferry 
fitting form that in a way is equivalent to the VFT form. 

Adopting the aforesaid TD-D correspondence established between the AGM and the 
VFT form given by Eq.(|^), and with the help of Eqs.(|^) and (|20|), we introduce the primary 
relaxation lengthscale through the relations 

(AG) ^oo ^ . , (AG) (VFT) ^-A-S* , . 

where the minimum cluster-size magnitude Za is anticipated by Eq.(P). 
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The CRR size-temperature behavior given by z!^^"^^ was experimentally studied for SCLs 
and for orientational-glass plastic crystals in Refs. [0,^]. These studies as well as the 
findings of Ref. |10], both based on the AGM, clearly indicate that Eg. (]23|) , and also Eq.(^), 
do good a job within the whole SCL ergodic temperature range, i.e., aXT >Tg (see e.g. Fig. 
4 in Ref. [^] and Fig. 3 in Ref. [1^). This permits one to present the experimental result for 
the excess isobaric heat capacity ^C^^^^ oc T~^, used above to evaluate the configurational 



entropy in Eq.(^), in the following equivalent form: 



(exp) 



^^^exp)^ for T>Tg. 



(24) 



With the help of Eq.(^ one arrives at a new useful relation for the observed model parameter 

introduced in Eq. (B|), namely 



A 5*00 (see Table I in Ref. 



A5oo = AC^^^^P) 



nig — nig 



(25) 



that though needs further experimental examination. 

Finally, one can see that the CRRs can be introduced by the SCL cluster statistics on 
the basis of Eq. ([18|) with the help of P^'^\z) = 6{z — z!^^^). This means that the CRRs 
are homogeneous clusters of characteristic size zt = z^r^'^'^ given in Eq.(4) with negligible 



heterogeneous and/or thermal size fluctuations, i.e., Az^ = Az 



(AG) 



0. 



C. Extension of the Vogel-Fulcher-Tamman Equation 



Within the cluster statistics approach the primary relaxation timescale is introduced by 
Eqs.([l5|),(p!^) and ([l8|) . A straightforward calculation using the standard integral 

„2 



exp ( — — — ax ] dx = V 7cb 



1 — erf ( a^/b 



yields 



1 + erf 



1 + AEr/kBT 
^y2AET/ET 



exp 



exp ba 



Ex I ^ /\Ex 



(26) 



2kBT 



(27) 



given in the pseudo-Arrhenius form with the cluster-energy characteristics, namely 



Et = ^^^TZT, AEr = ^1^1^^ = . (28) 

Zt i — io/J 



The standard error function that emerges in Eq. (|27|) is bounded by one, and ahhough it 
modifies the pre-exponential factor r^, it does not substantially modify the exponential 
temperature dependence of tt- Furthermore, the error-function term will be absorbed by 

Ta- 

We see that the VFT form of Eg. (|27D , when extended by cluster-size fluctuations, has 
the form of a high-T thermodynamic perturbation series. The first and second terms in 
the exponent describe the rigid solid clusters and their fluctuation corrections, respectively. 
The first term matches the rigid-cluster relaxation time given by the VFT equation and the 
AGM solution given by Eq.(^ and Eq.(^), respectively. The TD-D correspondence is now 
extended at a mesoscopic level that permits one to establish the physical meaning of the 
macroscopical phenomenological parameters for the strong-SCL states, namely 

j^iSL) ^ ^iSL) ^ A/VV^ .29) 

fee In 10 fcfilnlO ^ ^ 

under the additional condition A/i^ = A/i^"^*^) along with Eq.(|D. 

Unlike the case of the thermodynamic AGM, which gives a good description [1^ of all the 
SCL states (see Fig. 2) within the ergodic temperature domain, Tg < T < Tm (that means 
that C = = C(*^^) in Eq.(|l|)), the VFT form fails to describe weak SCL states 

if one uses the same critical temperature Tq. In order to incorporate the weak SCL states 
into the a-timescale parametrized by the observed MCT parameters {t^'^'^^ , Tc, and 7c) let 
us represent the exponential form ( p7D in the asymptotically equivalent algebraic form given 
by Eq.(|^. Using exp [c(l + x)] = exp [(1 + xY], one has the following high-T series with 
respect to the small fluctuation parameter AEx/ksT , namely 



(VFTE) _ 



AE' 

exp 1 + 



2kBT 



(30) 



AEt 1 / AEt ^ ^ 



(MCT) 



^ °°TN7c for T > T, (31) 
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denoted the VFT extended (VFTE) form. Accounting for only the two leading terms, one 
establishes the following relations, i.e., 7c = Ec/ksTc and = AEc/2kB with t^^'"'^^ = 
e'^'^Ta, that can also be represented with the help of Eg. (|28|) in terms of the observed MCT 
parameters, namely 



and 



A (VFTE) 

= To + . (33) 



Here a new model parameter = ^ AEt / Et = Azt/zt is introduced with the aim 
of establishing a scale for the cluster-size fluctuations above T^. Accounting for the fact 
that the size fluctuations are limited, at least to C,c < 1, with help of Eq.(|32D one arrives 
at a physical restriction 7c > 2 that is in accord with experimental data. As a result, the 
a-relaxation scale is now parametrized for the strong SCL states, namely 

4^^) = Ta exp i ^^^'f ] With = 1^ - 1, for T, < T < Tc, (34) 

with the observable combinations of cluster parameters Afigzi^^^ = DgksTQ (see Eq.(^)) 
and zT^^ = z^^'^'' (see also Eq.(||)). For the case of the weakly cooled states in moderately 
SCLs one has 



(ML) 
Tt =Ta exp 



for T,<T<Tm (35) 



ksToer y 7c kBT^er 

with A^cZit^^^ = IcksiTc — Tq) that follows from Eq.(|33D when z^^"^^^ = z^''^\ In 
the following sections we review experimental evidence for the proposed description of the 
primary relaxation through the timescale steepness and curvature. 



IV. STEEPNESS FUNCTION ANALYSIS 
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A. Qualitative Analysis 



We discuss two scenarios of relaxation of the SCLs observed in thermodynamic and 
dynamic experiments and shown in Fig. 2. The thermodynamic version is given by the AGM 
that describes the configurational entropy data within the whole SCL ergodic range (shown 
by solid squares in Fig. 2). Within the context of percolation treatment, the solid clusters are 
apparently the same for the strongly and moderately SCL states, i.e., zf^'' = Za^^^\= z^*^^) 
and thermal fluctuations are neglig ible (A4^^^ = 0). The AGM excess chemical potential 
is approximated by a constant given at Tg, i.e., A/i^^^) = Afig. 

To describe the primary relaxation timescale observed in the dynamic experiments, we 
start from the trial estimate given by the thermodynamic approximation. This is formally 
introduced by SCL-system clusters with zi^'^^^ = z^"^^ = Za, but with Az!^^^^ ^ 0, that 
corresponds to an extension of the AGM by cluster-size fluctuations, which strictly speaking 
are not observed in real thermodynamic experiments. Meanwhile the trial dynamic time- 
scale steepness function follows from Eqs. ([TO|) and (|35D , namely 



In 10 Et \ £t J \ £t J 



JJ{SCL) If ^ \ / 2 



j^iscL) ^ A/f^ Tg<T< T^, (36) 

where the auxiliary timescale function st is given in Eq. (^ . In order to employ the timescale 
parametrization in terms of the observed fragility m^, we apply Eq . (|36|) at T = Tg. Using the 
relation Za = kBTQ'jc^c/^fJ'c, one excludes the unknown parameter Za- By solving Eq.(^), 
with respect to the MCT parameter 7c, one obtains the trial estimate given by 

m* In 10 , , 

7c = / r 37 

A/ic eg \ A/lc Sg J 

that was analyzed in detail in Ref. Here we consider Eq.(P7|) as the quadratic equation 
with respect to the unknown factor A^g/Afic- This enables one to calculate the rebuilding 
energy per molecule near the crossover temperature, namely 



^ = /l + gll^lnlO-ll. (38) 
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The physical solutions of the quadratic equation (^TD must satisfy the requirement of ther- 
modynamic stability discussed above and given by the condition Afig > Afi^- Numerical 
solutions of Eq.(^) are shown in Table 2 for experimental data known for the charac- 
teristic temperatures and the slowing-down exponent 7^^^''^ given in Table 1. Predictions 
for the unknown 7^ are also made with the help of Eq.(P7D through the overall estimates 
AyUg'^'^^^/AyUc"^'"^^ = 1.2 aud A/i3^°'^/A/ii^°'^ = 1 adopted, respectively, for liquids and poly- 
mers. We see that no rebuilding chemical potential effects are found for the case of glass 
forming polymers. This can be understood as follows. In contrast to the case of molecular 
liquids, where application of the ergodic cluster description by the Gaussian distribution 
given in Eq.([T5D is justified both physically and experimentally, the monomeric-segment 



polymers require accounting of their fractal structure even for their ergodic states p9 

The dynamic relaxation scenario given by the VFT equation is realized through rigid 
clusters as well as by the VFT equation extended by the cluster-size fluctuations (VFTE 
equation) and applied, respectively, to the strong and weak SCL states. The relevant steep- 
ness function follows immediately from the corresponding relaxation times. Hence with help 
of Eq. ( P^ , one has 

m?^) = + -) with D, = (39) 

^ InlO^T St ^ knTo ^ ' 

for the case of the strongly SCLs. In the case of the moderately SCLs the steepness function 
follows from Eq. (p5D , namely 

^,M., = f 1 + L) I, + = (40) 

^ InlOET V ^t) V Ic^t) keTo ^ ' 

Qualitatively, the first-order derivative analyses of experimental data by Stickel et al., 
given in Refs. points to the existence of a small kink of the steepness function near the 
crossover temperature (see Fig.l) ||5^. Formally, this means that the limit of the timescale 
steepness function from below (T ^ Tc — and m^^^ — > mi^^^) differs with that from above 
(T — ^ Tc + and nij!^^'' mi^^^^ ), and therefore one should expect that m[^^^ > mi^^^K 
Application of these hmits to Eqs.(|39|) and (^OD using the relation Dc = ■Jc^c, that follows 
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from Eq.(^), results in Dg > SDc- These qualitative findings extended by the condition of 
the thermodynamic cluster stability can be summarized as 

^3Miy > 3 with — > 1. (41) 



The dynamic conditions given in Eq. (^T]) reflect the observed dynamical crossover from 
moderately to strong SCLs near Tc that are associated with the dynamic reconstruction of 
the cluster lengthscale. To provide deeper insight into the problem we will reproduce below 
the fitting derivative analyses by Stickel et al. within the VFT and VFTE equations. 

B. Numerical Analysis 

1. First- Order Derivative Timescale Steepness 

The first-order derivative analysis of the temperature behavior of the timescale was given 
in Refs. with the help of a first-order derivative function, namely 



IT 



(42) 



di/T J 

This function was estimated from r^*^^^^ and fitted by the VFT form, i.e. by F^j^^"^^ = 
(1 — Tq/T) / ^/B , and by other known fitting forms (see e.g. Fig.l). We reproduce this 
analysis in Figs. 5 and 6 in terms of the steepness function deduced from Eq . (|42|) , i.e., 
^(exp) _ i^j'(^p(^p)'-j2^ with F^y^'''* found by Stickel et al. by fitting with m^^^ and m^^^K 
The latter are given, respectively, in Eq.(^) and Eq. (^OD , with the strength indices Dg 
and Dc treated as adjustable parameters within the corresponding temperature regions, 
Tg <T < Tc and Tc <T < Tm, respectively. The results of fitting analyses are accumulated 
in Table 3. Unlike the case of Figs. 5, where the VFT temperature Tq is treated as the 
same for the both states, the SCLs analyzed in Fig. 6 are in conflict with the canonical VFT 
form for the strongly SCL states, but this is in accord with analysis given in Refs. [[7|,p|,p!0 



In this case the critical temperature Tq is not established and is therefore treated as an 
additional adjustable parameter for description of the moderately SCL states through the 
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VFTE equation. Thereby, we follow the idea that the Kauzmann temperature To might be 
established beyond the VFT equation. 

Our quantitative analysis of the strong-to-moderate-state crossover is in accord with 
the qualitative predictions given by Eq.(^). As seen from Table 3, the typical change in 
the ratio of strength indices Dg/D^. lies between 4 and 6. Exceptions should be given for 
the "irregular" case of salol where the Kauzmann temperature established by the VFTE 
equation (Tq = 2007^) is rather different from those established by other methods and given 
in Table 1. 

Adopting the trial estimate for the moderate-to-strong SCL-state crossover for the ex- 
cess chemical potential, given in Table 2, one arrives at the length-scale ratio estimate 

see that the observed crossover, associated with the cluster 
reconstruction, is not caused solely by decrease of cluster-size fluctuations and/or by growth 
of their chemical potentials. This crossover is accompanied by the lengthscale extension, 

i.e., .r^ > 



2. Second- Order Derivative Timescale Curvature 

The second-order timescale derivative-fitting analysis introduced by Stickel et al. can 
be treated as experimental observation of the critical temperature To through the VFT 
form. This analysis was done with the help of the second-order derivative function 
F2t{= 0(T)), introduced in Eq.(21) in Ref. as 

F.r^'-^(^^Y^T(^-'-^Y (43) 



dT \ dT^ J V d\nT ^ 

that results in ^2^^"^-* = — (T — To)/2 for the case tt = t^^'^^ given in Eqs.(|^) and @. 
The experimental data for rjf^^^ was used to estimate F^^^^ by Eq.(^) and fitted to the 
temperature-linear function F2^^'^\ As a result, the two distinct critical temperatures Tq 
and Tc were established for PC and 1-propanol through the VFT and HT-VFT forms (see 
Figs. 13 and 16 in Ref. M, respectively). In the case of salol and PDF, only the HT-VFT 



24 



form was derived in Figs. 8 and 11. Remarkably, in all cases the assumed kink in the 
function F^^^^ near Tc was transformed into the evident jump of the function ^2^^^"*. 

The curvature of the timescale kx, closely related to the second-order derivative function 
given in Eg. (^31) , can be introduced by the following relations, namely 

^ " dlnT ^ lnlO(rflnT)2 ^ ^ ~ F^' ^ ^ 

with k^^'^^ = (3T— To)/(T— To). In order to describe the weakly cooled states in moderately 
SCLs through the curvature k!^^^^ we substitute the relevant steepness function given in 
Eq.dD. This yields 

4t'^ = 7 ^ with = ^^fl (45) 



^ \lc+iDc ^ ET + l ) \ £T ) 

that leads to an expression for k!^^^^ with the help of Eq.(^4|). The typical SCL timescale cur- 
vatures are exemphfied by PC and salol given in. Fig. 7 and Fig. 8, respectively. We see that 
both cases are distinct and that the moderately SCL states are consistent with the cluster- 
size fluctuation description, based on the idea of the existence of a unique thermodynamic- 
instability temperature = Tq. They are described by the VFTE timescale equation (|35D 
with the same model parameters established in the steepness timescale analysis given in 
Fig. 5A and 6A, respectively, for PC and salol. Conversely, in the case of salol and PDE, our 
attempts failed |^ to the describe the strongly SCL states through the observed curvature 
with the help of various VFT-type fitting forms reviewed in Ref. including the proposal 
given in Refs. Accounting for the aforesaid TD-D correspondence, we speculate that 

the failure of the standard VFT equation in the salol-type SCLs (that should be extended by 
the case of BMPC (see Fig.6C)), is due to the inapplicability of the cluster description given 
by the thermodynamic model of Adam and Gibbs in Ref. |l^. Indeed, as shown above, 
the main ACM constraint for the excess chemical potential, i.e., AfiT = ^fJ'c, works well 
above the SL-ML crossover temperature, but below it needs to be reconsidered for the 
salol-type SCLs. 
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V. DISCUSSION AND CONCLUSIONS 



The investigation of the temperature behavior of the primary relaxation timescale in 
SCLs plays a fundamental role in understanding the glass transition problem in general, 
as widely discussed in various, microscopically distinct, glass forming non-solid materials. 
In spite of the fact that the nature of the structural glass transformations is essentially 
dynamical, the process of solidification of rearranging regions, while avoiding crystallization 
of a given SCL as a whole system, exhibits many features of the common liquid-to-solid or 
disorder-to-order true thermodynamic phase transitions. 

The "pure" dynamic treatment of the glass-transition problem, given by the simplified 
MCT version, and the "pure" thermodynamic treatment given by the simplified AGM, are 
limited in the description of the primary relaxation timescale temperature behavior. This 
limitation was revealed by the temperature derivative analyses given by Stickel et al., as 
shown in Fig.l, that clearly establish that the fitting forms for the known t!^^^^^ and t!^^^ 
(given, respectively, in Eqs.(|^) and (0)), fit the data only for a restricted temperature 
window. 



The same kind of restriction applies to the analysis of the thermodynamic models |12|JT3 
mentioned in the Introduction, which offer a similar overall r-p description within the whole 
range of existence of the ergodic SCL states. Although these models seem to provide a good 
fit to 4^=^^^ (by using only one adjustable parameter, as the fragility nig in the case of Ref. 

or by using more than three adjustable model parameters, as in the case of Ref. [|Tl[] ), 
they need to be analyzed by considering the derivative functions of the dynamical variables 
over the full temperature range. 

As stressed by Angell et al. in Ref. 0], the popularity of the thermodynamic AGM 
among molecular-liquid phenomenologists might be compared with that of exploring the 
phenomenological VFT form by glass-material experimentalists. The current study develops 
a statistical approach for the description of CRRs, with their isobaric-isothermal mechanism 



of solidification that were introduced a long time ago by Adam and Gibbs in Ref. |]T3|. 
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The application of the second-order derivative analysis for the case of salol, revealed 
that the widely used standard VFT form given in Eg. ([T6|) has no observation window. 
This observation ended the unreserved domination of the phenomenological VFT form for 
fitting all dynamical data of all SCLs, and divided the SCLs in two groups, "regular" and 
"irregular", that, respectively, obey or do not obey the VFT equation. Meanwhile, it seems 
that all SCLs exhibit a jump in the timescale curvature near Tc (see e.g. Figs. 7 and 8) 
related to the corresponding kink observed in the temperature behavior of the timescale- 
steepness that separates the glass forming liquids into two dynamically distinct states (see 
e.g. Figs. 5 and 6). Additionally, both the "regular" and "irregular" SCLs, therefore manifest 
the absolute thermodynamic instability at the Kauzmann temperature T^, whose physical 
meaning is extended by the equation = Tq. 

We have discussed the SCLs in ergodic states in terms of isolated finite-size solid-like 
clusters observed within the window Tg <T < Tm in the dynamic experiments through the 
primary relaxation timescale data rjf^''^ The cluster dynamics, introduced by the Gaussian 
size-distribution cluster statistics and exponential- type size-relaxation, described by Eqs.(|T3p 
and (0), respectively, is driven by the mesoscopic entropy and ensured by the Onsager 
thermodynamic principle. 

Below the crossover temperature Tc, the CRRs emerges as structurally homogeneous 
clusters, where thermodynamic stability is guaranteed by the growth of the excess chemical 
potential A/i^ (> A/ic). This energetic change, associated with the cluster rebuilding, is 
justified experimentally through the analysis given in Table 2A. In turn, the homogeneous 
clusters are associated with the strong SCL states of the MCT, and/or with percolation 
theory solid-like clusters described by the VFT equation (with the activation energy Ej, = 
A/igZ^^^^/T and the characteristic size zf^^^ = zi^^'^ /{I — Tq/T), both are given in Eq.(|3iD). 

(SL) 

The minimum cluster size Za is related to the lattice constant of the underlying percolation 

(SL) ( (SL)\ ^^'^ 

lattice, i.e., Q = ( voZa ) , where vq is the molecule volume. 

The experimental observation |TD| of the strong-SCL ergodic states, described simulta- 



neously by the thermodynamic ACM and the dynamic VFT equations led to formulation 
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of the thermodynamic-dynamic glass-transition correspondence principle. This principle is 
based on the physically intelligible thermodynamic-dynamic equations introduced for the 
critical Tq and glass transition Tg temperatures, as well as on the ergodic hypothesis that 
gives rise to the temporal and spatial scale minima introduced in Eq.(|^). These scales are 
common for both dynamic and thermodynamic experiments within the temperature domain 
Tg < T < Tc- Upon cooling near and below T^, the cluster size z.^^^^ and its heteroge- 

ISL) 

neous fluctuations Azj, grow and the SL clusters show an ergodic instability caused by 
their maximum fluctuations |2^. The compact "ergodic" solid-like SL clusters are there- 



fore transformed into fractal-type clusters of the holey-like (approximated by the ideal gas 
p7| , |29[| ) structure that is accompaniment by expansion of the relevant percolation lattice, 
with Cfe > Ca • Upon heating, the CRRs of the AGM disappear in the dynamic experi- 
ments above but remain stable in calorimetric thermally equilibrated experiments (see 
Fig. 2). Moreover, as follows from the first and second order T-derivative analyses elaborated 
in Refs. [[7|,p|, p^p5[ , the SL-cluster description given with the help of the VFT-type forms 



is not useful for "irregular" SCLs such as salol, PDE, and BMPC (see also steepness and 
curvature analyses in Figs.6A and 8, Fig.6B, and Fig.6C, respectively). 

Unlike the SL states, the structurally heterogeneous clusters, denominated as ML clusters, 
exhibit stability in all studied SCLs. They are associated with excitations of the weak 
supercooled states, characteristic of the moderately SCLs within the domain Tc < T < 
and introduced by the VFTE equation with the correlation length zf^^^ = Za^^^^/{1 — 



(ML) 
T 



Tq/T) described in Eq. (^8]). In turn, the heterogeneous cluster-size fluctuations, Az, 
(2/7c)^/^ z^^\ are introduced with the help of the model parameter given in Eq.(^). 
Due to the common physical restriction imposed on fluctuations, which preserve the ML- 
cluster structure and given by Azi^^^'^ / z^^^'' < 1, one naturally arrives at the constraint for 
the timescale slowing-down exponent, i.e., 7c > 2. 

The minimum size of the ML clusters derived from the dynamical experimental data (see 
Table 3) is related to the instability of heterogeneous clusters below the crossover tempera- 
ture Tc that results in an expansion of the underlying percolation lattice, i.e., C,^^^ > 
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We see that within the context of percolation theory, the discussed cluster reconstructions 
near Tg and Tc are geometrically similar and therefore the crossover temperature for the 
moderately SCL states plays a role similar to the "crossover" temperature Tg for the strong 
SCL states. This clarifies the effectiveness of the phenomenological HT-VFT fitting form 
introduced on the basis of the VFT equation with relevant substitution of the characteristic 
temperatures discussed in Eq.(|l^). 

It has been justified experimentally (see analyses in Figs. 5-8 ) that the size fluctuations 
of the ML clusters, as "non-rigid" solid-like precursors of the SL clusters, grow upon the 
decrease of the temperature. This means that the ML clusters are characterized more 
as "fluid-like" rather than "solid-like" and hence might be introduced on the basis of the 
thermodynamic heterophase-fiuctuation microscopic model proposed by Fischer and Bakai 
in Ref. According to this model, the ML clusters may be associated with isolated liquid- 



like droplets that show small heterophase fluctuations. Upon cooling, the concentration of 
the solid fraction grows and the liquid droplets exhibit a crossover from isolated to non- 
isolated long-range density fluctuations |Q, related to Fischer clusters that eventually do 
not contribute to the primary relaxation. As seen from Fig. 5 in Ref. P2[, for the case of 



salol, the droplet percolation crossover temperature is close to i.e. to that predicted by 
the MCT. The question therefore arises as to whether the interpercolation droplet states. 



given by the Gibbs potential in Eq.(2.1) in Ref. ||22|, enable one to describe the drastic 
temperature behavior observed for the timescale curvature in the window Tg < T < 
Tc (see Fig. 8). Furthermore, one can speculate that the thermodynamic-dynamic glass- 
transformation correspondence above T^, can be established on the basis of the physical 
correspondence between the thermodynamic heterophase-fiuctuation droplet model p2|] and 
the dynamic VFTE equation given for t!^'^^^ in Eq. (pSj) . 



T 

One can expect that application of the derivative analysis by fitting the timescale steep- 
ness and curvature, in addition to the common t^^^^ fitting, might serve as a reliable in- 
strument to give new insights into important features of the underlying physical processes. 
This is related to the weak-to-strong SCL-state crossover based on the existence of the 
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observation windows and the steepness kink discussed above in Eqs.(|^) and (^T]), respec- 
tively. We exemplify our statement by the Schulz model (SM) |12[ that employs the concept 
of Goldstein's energy landscape for SCLs. The overall ergodic-range timescale |12[ can 
be represented here by the activation energy E^^^^ = Em{Tm/T)"'^ where Em and 7^ are 



m 



thermodynamic parameters (with 7^ = 1 + 2Ac*-'^*^Y3^b , given in Eq.(26) in Ref. Wl 
terms of the excess liquid-over-HT-solid isobaric specific heat capacity Acf^^^). Ignoring the 
weakly T-dependent pre-exponential factor, one can estimate the steepness function defined 
in Eq. ([lO|) at Tg, i.e., the fragility nig. Similar to the case of the VFT form, discussed on the 
basis of Eq.(|^), one can solve the system of the two equations for rrig and Tg with respect to 
Em and ■jg. This process results in the following expression for the slowing-down exponent, 
1g ~ ''^g/^l ) cind offers a new equation for the Schulz model excess heat capacity, namely 

We have thereby demonstrated, along with the findings given in Eq.(0) and analyzed in 
Fig. 4 (for the self-consistent characteristic temperatures) and in Eq.(P^) (for the observed 
excess entropy A^oo), how the steepness timescale analysis can be used to verify the self- 
consistency of the model predictions. Both model Eqs.(^) and (^) need further testing. 

We also speculate on the existence of some alternative mathematical forms for the 
timescale tt, whose cluster-size singular behavior near Tq is sufficiently affected by the 
excess chemical potential A/it- Indeed, in the current study, the macroscopic primary 
relaxation timescale is introduced by the mesoscopic solid-like clusters with rearranging 
frequency u!^^\z) = t^^{z) = raexp(— AyUy z/ksT) (see Eq.(p!7D). An alternative descrip- 
tion of the random mesoscopic rearranging regions was proposed by Odagaki in Ref. 



on the basis of the AGM treated within the controUed-diffusion approach. Formally, the 
rearranging regions in SCLs were introduced as random domains defined by the random 
excess chemical potential AyU and by the characteristic size z^'^'^ that follows from Eq.(pTD 
given by the AGM. The domain rearranging frequency was therefore treated as a random 
value a;^°™'^(A/i) = t^^ exp(— AyU z!^^^ /ksT). In fact, the most general mathematical form 
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is given by ujt = Tr^i^l^^ z) = exp(— A/i z/ksT) where both A/i and z are random 
variables. This type of extended generahzed treatment might be considered as a possible 
candidate to solve the problem of the primary relaxation observed in the salol-type SCLs. 

Finally, we have demonstrated how the primary relaxation scale can be described by the 
two VFT and VFTE equations and parameterized by a finite number of observables: three 
principal characteristic temperatures Tq, Tg and Tc; the fragility m^; the slowing-down expo- 
nent 7c; and the two strength indices D^^^^ = Ajjgzi^^^ /kBTo and D^^^^^ = Ajj^^zi^'^^^ /ksTo 
given in Eqs. (|39|) and (0), respectively. Physically, this description is justified by the sta- 
bilization of the strong and weak SCL states, realized through stabilization of the temporal 
and spatial scales, and related to the minimum time and size of the rearranging regions. This 
proposed approach is elaborated with the help of a set of "universal parameters" m* = 16±2 
and m* = 7±1, that seems to be extended by the new ratio zi^^'^ / zi^^^'^ = 4.2±0.5. Further- 
more, one can show through application of the TD-D correspondence principle, that CRRs 
are observed with the minimum configurational entropy, i.e., Asg^^^ = s* = ln(2 ± 0.25), 



as predicted by Adam and Gibbs in their pioneer work in Ref. The new result will be 
discussed elsewhere. 
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FIGURE CAPTIONS 



Fig. 1. First-order derivative function of the timescale, 10^ ((ilog^Q(rT)/(i(l/T))^^^^ ver- 
sus temperature for propylene carbonate. Circles: relevant dielectric relaxation and dc 
conductivity experimental data r^^^^^ from Fig. 12 in Ref. Lines: fitting forms proposed 
by the mode coupling theory (MCT), Adam-Gibbs model (AGM) and phenomenological 
VFT (To = 132i^, B = 389K) and the HT-VFT {T^"^^ = 153K, S^^^) = 158K) forms. 
One can see that the AGM and VFT forms fit the strongly SCL (SL) states, shown between 
Tg = 157K and = 17QK, and the HT-VFT form fits the moderately SCL (ML) states 
and even normal liquid (NL) states above the melting point = 225K . Ty is the lower 
temperature of the MCT form. 

Fig. 2. Typical temperature behavior of the configurational entropy of supercooled liq- 
uids. Schematic representation of the thermodynamic data [|T^ (solid squares), reestimated 
dynamic dielectric relaxation data (open circles, see Eq.(8) in Ref. [|Tn|), and the mechanical 
stress relaxation data (solid circles, see Fig. 10 in Ref. [^) data. The solid hne corresponds 
to Eq. (|^), which is common to the AGM and VFT forms for strongly SCL states, and fits 
the thermally equilibrated states within the whole ergodic region. Notations correspond to 
Fig.l. 

Fig. 3. Observation of the universal parameter m*: scaled characteristic temperatures 
against scaled fragility for glass forming liquids (A) and polymers (B). Circles (A) and 
diamonds (B) correspond to experimental data, see Tables lA and IB, respectively. Solid 
line obeys Eq.(§) with m* = 16. 

Fig. 4. Experimental justification of the characteristic-temperature equations for super- 
cooled liquids (A) and polymers (B). The symbols denote experimental data from Table 1. 
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Solid and broken lines are solutions of the CTEs, ([I3D with fixed ratio m*/m*, and (^4|), 
respectively. Dotted lines correspond to the upper and lower limits of Eq.(p!3|) subject to the 
condition m*/m* = 2. 

Fig. 5. Steepness-temperature analysis of the primary timescale in supercooled liquids 
propylene carbonate (A), ortho-terphenyl (B) and 1-propanol (C) where the VFT equation 
(P) is experimentally justified. Symbols correspond to the dynamical data (dielectric relax- 
ation, dc conductivity, and viscosity) by Stickel et al. reestimated with the help of Eq.(^2|) 
from Fig.l2 in Ref. P], from Fig. 6 in Ref. and from Fig. 15 in Ref. |^ for the cases A, 
B and C, respectively. Lines are given by the VFT and VFTE equations using the experi- 
mental characteristic temperatures from Table 1 and with the adjustable parameters given 
in Eqs.(|39D and (j^, respectively. Notations correspond to Fig. 1. 

Fig. 6. Steepness-temperature analysis of the primary timescale in supercooled hquids sa- 
lol (A), phenolphthalein-dimethyl-ether (B) and bis-methoxy-phenyl-cyclohexane (C) where 
the applicability of the VFT equation (||) is questioned. Symbols correspond to the dy- 
namical data given by Stickel et al. and reestimated with the help of Eq.(^2|) from Fig. 4 
in Ref. 0, from Fig. 10 in Ref. and from Fig. 7 in Ref. |^ for the cases A, B and C, 
respectively. Lines and fitting parameters correspond to Fig. 5. The VFT fit is questionable 
at low temperatures. 

Fig. 7. Curvature-temperature analysis for the primary relaxation in propylene carbon- 
ate. Symbols correspond to experimental data reestimated using Eqs.(|^) and (|45|) with the 
help of Fig. 13 in Ref. 0. Lines are given by the VFT equation {k!^^'^^ ) for the strongly 
SCL state (SL), HT-VFTE (4^^^^^^^ ) and VFTE {k!^^^^^ ) equations for the moderately 
SCL state (ML), all discussed in the text, with the same adjustable parameters as found in 
Fig. 5A. 
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Fig. 8. Curvature-temperature analysis for the primary relaxation in salol. Symbols 
correspond to experimental data reestimated from Fig.8 in Ref. [0. Lines represent equations 
given in the text with the same adjustable parameters as found in Fig. 6 A. 
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